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ABSTRACT 

'^q' . We give the modern derivation of Deser's formula using analytic continuation of the 

Q.^ ' scattering amplitude as a function of momentum. The electromagnetic corrections to the 

I pionium lifetime are given as evaluated in a potential model. 
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1. Introduction 



The Deser formula, first given in Ref . [1] , relates the low energy scattering properties 
of strongly interacting charged particles to the hadronic properties of the atoms built from 
! them. The main examples are: 

^ I n~p elastic and inelastic scattering and vr-mesic hydrogen (vr~p): 

l> , 

O ■ TT p — »• TT p, 7r°n or with (vr p) Tc^n or jn 



71 7r+ elastic and inelastic scattering and pionium (vr 7r+): 

— > 7r~7r^, tt^tt" or 77 with (n'li^) n^n^ or 77 



Oh! 

r-| ! To be more specific from the start we concentrate on the tt^tt'^ system. Nearly all 

information comes from the s-wave state, so we treat only / = and do not indicate this 



by a specific label. The 77 channel does not play a decisive role, so we forget it for the 
moment. In principle it could be included in the analysis. 
5^ ! If we compare physical 7r~7r"'" scattering with pure Coulomb scattering of two op- 

positely charged particles we measure certain scattering lengths; if we compare the bound 
state (7r~7r+) with the pure Coulomb bound state of such an atom we measure an energy 
shift and a finite lifetime of the ground state. 

The aim of the Deser formula then is to relate the scattering lengths to the shift 
and decay width of the ground state of the atom. For the two channel scattering system 
we have two scattering lengths acc and aoc : 



_ yTT 7r+ : ttcc 

We will give the precise definitions of these quantities later. It will turn out that 
the shift of the ground state is dominated by a term proportional to and that the 
width for the decay 
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is dominated by a term proportional to |aoc|^. To understand the main ideas for the 
derivation of this result it is convenient to treat a simplified system first: We assume 
that the 7r° does not exist. In this case only elastic scattering would be possible above 
threshold and pionium would be stable; this one channel problem will be treated in Sect. 2. 
In Sect. 3. we will then treat the realistic two channel case and in Sect. 4. the problem of 
isospin conservation and electromagnetic corrections to this symmetry will be dealt with. 

2. One channel case 

As already described we assume that 7r° does not exist. We only deal with elastic 
7r~7r+ scattering and the bound state, which in this case could not decay for energetic 
reasons. But due to the hadronic interaction the energy of the bound state will be shifted 
compared with the energy of the pure Coulomb ground state. So we face the problem of 
connecting the energy shift with the scattering length for 7r~7r"'" elastic scattering. The 
problem has been solved in a semiquantitative way in Ref. [1]. A systematic method 
which also allows the inclusion of corrections has been worked out on the basis of analytic 
continuation of the scattering amplitude in Ref. [2,3]. To this end we have to distinguish 
between 3 different scattering amplitudes: f^^^ {H for hadronic), f^'"''^ {Ch for pure 
Coulomb) and f^^^ (T for total). In the following q is the magnitude of the momentum 
of either particle in the cm-frame; rric is the mass of the tt^. 

We will now treat the different scattering amplitudes in a way which later allows 
an easy generalisation to the two channel case. 

-purely hadronic scattering 



fW-i^q)=K^^^-\q')-tq. (1) 

K^^") is a real function of g^, analytic near = 0. 

One can expand K^^^'^^q"^) in the following way: 

K^^^-\q')^a(^^-' + lr^f^}q' + 0(q% 

a^^^ is the hadronic scattering length, i.e. the scattering length with the Coulomb inter- 
action neglected; r^J is the corresponding effective range. 



2// 

From Eq.(l) we have: 



,(«)(„) - ^'"'fa^ (2) 

A pole of f^^^ on the positive imaginary axis of the complex q - plane corresponds to a 
hadronically bound state. The only possibility for such a pole near q = can come from a 
zero of the denominator in Eq.(2). In our case there is no such hadronically bound state. 

-pure Coulomb scattering 



The corresponding scattering amplitude (partial wave amplitude for Z = 0) is known 
explicitly: 

•' 2iq ' 
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with 

r(l-ir^)' 



am 

V = 



2q qas 

Here a is the fine structure constant and is the Bohr radius of pionium. Prom the well 
known analytic properties of the F-function one finds the poles of f'^'^^^ (q) on the positive 
imaginary axis of the complex q - plane at the positions 

i 1 

qn^ , n = l,2,3, ... (3) 

n qb 

They correspond to the Coulomb bound states for I — 0. 
-Mixed scattering 

We assume that there is a short range hadronic interaction as well as the long 
range Coulomb interaction, which are additive in the Hamiltonian. The resulting total 
scattering amplitude then is 

y(r) ^ y(Cb) _^ g2iao 

with 

f-\q) = C^\ri){K-'{q') - 2qnh{ri)) - iq, (4) 

where 

Co' = e-'^T(l + i77)r(l-i77). 

h{ri) also is an explicitly known function, related to the digamma function. It is 
well known that K{q^) has the same analytic properties as K'^^\q^) near q^ = 0: 

K~\q')=a-' + \r,ffq' + 0{q% 

Here a is the scattering length in the presence of the Coulomb potential and Tg// 
the corresponding effective range. We then have 

^ ^2iq ^ 1 - K{q^){2qr^h{r^)+iqCl{r^)y 2iq 

We are now looking for the poles of /^^^ as a function of q. e^"^° has (as we know already) 
poles at qn (see Eq.(3)) for the pure Coulomb states. A careful analysis shows that the 
curly bracket vanishes at g„, so the product does not have a pole. The only poles of f^'^^ 
therefore can come from the solutions of (note that rj depends on q) 

1 - K{q'^){2qrjh{r]) + iqC^{v)) = or equivalently r\q) = 0. (5) 

These lie on the positive imaginary axis, near the pure Coulomb poles; the differ- 
ence in the position gives the energy shift AE due to the hadronic interaction. The result 
for the Is state is: 
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AE = —{1 + CI— + ...)Ei, (6) 

where Ei is the pure Coulomb ground state energy corresponding to qi = i/aB- 

If one replaces a by a^^^ and neglects the second term in the parenthesis of Eq.(6) 
one finds (in our notation) the result of the semiquantitative treatment in Ref.[l]. The 
second term in the parenthesis is (since a/a^ <^ 1 ) numerically negligible compared to 
the accuracy expected from experiments. It contains the higher order hadronic effets on 
the bound state. Furthermore ci depends on Vcff] this takes into account the subthreshold 
effect, which means that K{q^) has to be continued to negative values of q^. 

The difference between a and a*^^^ is well known from the comparison of pp- 
scattering with np-scattering in the I — S — T— 1 state of the two nucleon 
system. In this case a^^^ (measured from np scattering (no Coulomb force!)) is approxi- 
mately 2.5a (a being measured from pp scattering). Thus the difference exists and can in 
principle be measured. In the two nucleon system it is enhanced by a so called "virtual" 
bound state just above threshold. In the 7r~7r"'" system it is fortunately not as big as that 
and can be treated by perturbation theory. Since it is an electromagnetic effect we expect 
a difference of a few percent. 



3. Two channel case 



In reality we have to treat the coupled two channel case (still neglecting the 77- 
channel) : 



Let q be the magnitude of the momentum of each incident particle in the center of mo- 
mentum frame, go(Q') the magnitude of the momentum of each outgoing 7r° in that frame 
(mo is the mass of tt"): 



2 I 2 2 I 2 
qQ + mQ = q + 



Again we treat only s-wave scattering, which is justified at very small energies. The 
scattering amplitudes /^'^ now become a symmetric 2x2 matrix: 



F(-)(g) 



Fi-hq) = Fi^{q) F^iq) 



The matrix index c refers to the 7r~7r+ channel and refers to the 7r°7r° channel. We will 
treat this problem first by neglecting the Coulomb interaction in the scattering problem 
(i.e. above threshold). But we keep the option of treating rric 7^ mg. We consider the two 
cases: 

I : mo = rric ; II : mo < rric 

I: Since besides the Coulomb interaction we have also neglected the mass difference we 
are in the hadronic situation and indicate this by a superscript (H) in the scattering 
amplitude and scattering cross sections. We are dealing with coincident thresholds and 
therefore F'^^^O) — a^^\ with real matrix elements a^\ a^^^ = o^oP ■, <^m^ as in the 
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single channel system. The hmiting values of the tt 7r+ elastic and charge exchange cross 
sections are: 

a) lima(f)(g) = 47r|Fif )(0)r = 47r|aWr 

b) limaSf (5) = 47r|F(f (O)l^ = 47r|4f ^ 

II: Here the mass difference is taken into account and we indicate this by a superscript 
(A). The matrix F(^)(0) now has complex matrix elements. The limiting behaviour of 
the cross sections is: 

o) lim<T(fM = 4x|Fif>(0)|^ 

Going over from la to the more realistic situation Ila is not a dramatic change. We can 
expect that F^^\0) and a^^^ differ by only a few percent. But going over from 16 to 116 
changes the low energy behaviour of a^f^ (g) drastically, because we are now dealing with 
a real exothermic reaction. Nevertheless we can expect that Fq^\o) and aQ^"* also differ 
by only a few percent. Going from I to II obviously breaks a possible isospin invariance 
of the hadronic interaction. 

We now treat case II more systematically. As in the single channel case we have 

F(^)-^(g) = i^(^)-i(a-^g(g), (7) 

where 

- ( J,) ) ■ - ( Jo) ) 

From unitarity one can show that K^^\q^) is a real symmetric 2x2 matrix and that in 
analogy with the single channel case 

For our present pedagogical purpose we only need Eq.(7) at g = 0: 
Solving for the matrix elements of ^'•^•'(O) we get: 



F4^)(0) = 4^)+.go(0) 



1 - i?o(0)4o ^ 



(A) 



For(0) = 



l-^go(0)4?^' 

In the limit mo = rric we have qo{0) = and F^^\0) = yl(^)(= a^^)). 

We see that, in contrast to the one channel case, the threshold elastic scattering 
amplitude ^^'^^(O) now has an imaginary part. Replacing a in Eq.(6) by F^^\0), we have 
(neglecting terms of higher order in the mass difference and 0{a/aB))'- 



AE(-) = l^E, + z4,o(0)(4tO'- (8) 
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Here the real part of AE^^^ gives the energy shift of the ground state as found before. 
The imaginary part of AE^^^ is connected to the hne width F and the hfe time r by 

/mAE(^) = -ir(^); ^(a) ^ r^^)-^ 
2 

We therefore have 

r(^)- = -8,o(0)(4^))'^. (9) 

This phenomenological derivation of Deser's formula for the two channel case has been 
given in Ref.[4]. 

The final systematic derivation of Deser's formula has been achieved in Ref.[5]. It 

includes the Coulomb interaction also above threshold in addition to the mass differences 
and makes use of the method of analytic continuation in the two channel case. In analogy 
with Eq.(4) a matrix F{q) can be defined by 

F~^{q) — K~^{(f') — iQ{q) + Coulomb terms in the charged channel, (10) 

where K{q'^) again is a real symmetric matrix analytic near q'^ — 0. We define K{0) — 
A. We do not give the Coulomb terms explicitly nor do we comment on the analytic 
continuation as we did in Sect. 2. The result is that we have (analogously to Eq.(5)) to look 
for the zeros of detF^^{q). These zeros are off the positive imaginary g-axis, corresponding 
to complex energy shifts of the Coulomb poles as in Eq.(8). The interpretation is the same 
as given after Eq.(8). We only give the full result for the lifetime r of the ground state: 



1 8 go(^c)' ^ , ^ 



El [l + O [ — ^ ) + subthreshold corrections \ (11) 



r aBl + qliAoaY \ KaBJ J 

Here go is the magnitude of the momentum of each 7r° in the rest frame of the pionium 
atom. Neglecting the numerically unimportant subthreshold corrections as well as the 
terms of 0(/loc/a_B) and O {(Iq^AqqY) we get nearly back to Eq.(9); but note the difference 
between % and go(0) and the different definitions of A^^^ and Aqc- 

We emphasize very strongly that the Deser formula connects measurable quantities 
with measurable quantities: The matrix elements of A can "in principle" be measured 
in low energy tt~tt^ elastic and inelastic scattering and tt^tt" elastic scattering. The 
expressions for the cross sections in terms of A^c, ^oo; including mass differences and 
the Coulomb effects near threshold are given in Ref . [6] for the 7r~p - system and can easily 
be translated to the t:~t:'^ - system. 

4. Isospin invariance and electromagnetic corrections to the pionium life time 

Up to now we have not made any use of a specific model for the hadronic interaction 
nor did we assume a symmetry property like isospin invariance. 

From now on we use the well known result of chiral perturbation theory that the 
hadronic interaction in the tttt - system is isospin invariant. Whilst the proper under- 
standing of radiative corrections to this result needs a calculation in the framework of 
chiral perturbation theory Ref. [7], an estimate of their order of magnitude can certainly 
be achieved in a simple phenomenological model. Such a model has been developed in 
Ref. [8,9] for the 7r~7r+ - system. It assumes that the bulk of the electromagnetic effects 
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has its origin in the mass difference between tt" and and in the Coulomb interaction 
between vr" and 7?+. This method has been applied in the past to the pp - system and to 
the 7i~p - system (Ref.[5,6]). 

We will concentrate here on the application to the lifetime of pionium; so we are 
interested in the quantity 

Aoc-aS\ (12) 

with 

aiT-%a^-a'), (13) 

where and a° are the hadronic values of the s-wave scattering lengths for total isospin 
2 and respectively. These two quantities come from chiral perturbation theory without 
electromagnetic interactions and can be evaluated either for the reference mass mo or for 
the reference mass rric. 

It is well known from chiral perturbation theory that the hadronic mass of the pions 
is practically equal to mo and that the difference rric — mo is of electromagnetic origin. 

( 

Therefore taking mo as the reference mass for evaluating does result practically in a 
purely hadronic value. On the other hand, taking rric as the reference mass results in an 
isospin invariant, but not purely hadronic, value for aQ^\ Neverthless, as in most of the 
chiral perturbation theory work, we take m^ as the reference mass for the following. 

We emphasize also that we calculate the electromagnetic corrections to the hadronic 
amphtudes and that the values of the corrections are hardly affected at all by the uncer- 
tainty in the starting values of and a° (one loop result, two loop result, generalized 
chiral perturbation theory, ...). 

We now describe how the potential model is used to estimate the electromagnetic 
corrections to the scattering lengths for physical processes. Of special interest is the 
difference in Eq.(12). Our input is the isospin invariant hadronic scattering phase shifts 
S^{q) and 5°(?) given by chiral perturbation theory with the reference mass rric- The 
isospin invariant scattering lengths are then given by 

= lim (5^ (g)/g) T = 0, 2 

Taking as an example the results of Ref.[10] we get: 

a° = 0.2883 fm, = -0.0617 fm 

and from Eq.(12) 

4f = -0.1650 fm. 

In each isospin channel T we then construct a potential which inserted in a relativisti- 
cally modified Schrodinger equation with reduced mass |mc reproduces these phase shifts 
and scattering lengths. 

The two potentials V'^ are then used in a relativistically modified coupled two 
channel Schrodinger equation in which the Coulomb potential is added in the charged 
channel and the physical masses rric, "^o are used. Solving this two channel scattering 
problem we get the matrix K{q'^) (see Eq.(lO)) and A — K{0). Thus we have the matrix 
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element Aq^ which dominates the Deser formula Eq.(ll) for the lifetime of pionium. The 
result of this calculation is 

Aqc - 4? = 0.005(l)fm (14) 

The 20% error on the electromagnetic correction in Eq.(14) reflects the fact that the 
construction of an energy independent V'^{r) from S^{q) is not unique. We have tested 
several forms for such a potential which reproduce the phases well up to 500 MeV cm 
energy. We have also tested energy dependent potentials V"^. In all these cases, the 
difference in Eq.(14) lies within the error given. 

Numerically we arrive at the following relation connecting the hadronic quantity 
Qq^^ given in Eq.(13) and the lifetime of pionium r : 

0.0951 

r(fs) = 2 

(4f(fm) + 0.005) 

{a^c^ evaluated using the reference mass rric)- 

We have emphasised explicitly once more that rric has been taken as the reference 
mass. If niQ is used instead, the potentials V'^ have to be determined from a Schrodinger 
equation with a reduced mass |mo and the electromagnetic correction has a different 
value. This point of view was taken in Ref.[8], where the corresponding results can be 
found. 

We also emphasise that the influence of the 77 channel on the electromagnetic 
corrections to the scattering lengths can be incorporated by doing a three channel analysis. 
The analogous treatment of the 772 -channel in the 7r~p -system has been performed in 
Ref. [5]. Numerically this additional correction is negligible. 

We thank W.S. Woolcock for a very careful reading of the manuscript. 
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